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Previous speakers have discussed how to change the analytical framulation of various prob- 
lems in order to introduce stabilizing effects or to imivove the accuracy and thereby in- 
crease the efficiency of numerical integration. The present topic is concerned with how to 
improve the process of numerical int^ration itself, in m'der to increase efficioicy in the 
development of accurate ephemerides for earth satellites. In particular, the subject of this 
discussion is a new class of linear multistep methods for the rumerical integration of ordin- 
ary differential equations. These methods are distinguished from the classical methods in 
that they permit the solution to be corrected at cotain “back” points. That is, in the case 
of satellite computations, the solution is corrected at certain points in the past as the 
int^ration advances in time. Algmithms have beoi developed for the solution of both 
first- and second-order differential equations, althou^ only the second-order case is con- 
sidered here. 

There are two reasons for correcting the solution at back points: First, when a polynomial 
obtained from interpolating evenly spaced data is used for approximating a function at a 
point, the coefficient of the errm- term is smaller when the point is nearer the middle of the 
range of data. Therefore, by performing the last correction at an internal point of the grid, 
a more accurate Isolution (that is, smaller truncation error for a given order) is obtained. 

The second (and not so intuitive) reason is that the introduction of back corrections induces 
numerical stability. It is well known that the general stability boundaries for the predictor- 
corrector methods of the Stormer-Cowell type decrease geometrically with increasing order. 
That is, when the higher order Stormer-Cowell methods are used on a practical problem, 
the step size is severely constrained by stability considerations. This is not the case with the 
methods using back corrections. In fact, in some cases, these methods possess general stab- 
ility boundaries which are 30 to 40 times larger than those of the Stmmer-Cowell methods of 
the same order. Moreover, the stability regions of these methods do not exhibit geometrical 
decay but remain relatively unchanged up to the eighteenth or nineteenth order. 
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The back<orrection methods are given by: 

k 

■•‘2) ^ -(m + 1) ^ or. f . 

^♦1 ^ ' n«m ' ' n-iii>>l j n-j 

j=0 

and 

k 

^+i-« = (n» + 2 - 8) x„.^ -(m + 1 - 8) Pjb 

j=o 

for 8 = 0, 1 , .... m. The position vectcnr of the satellite at time t^ is denoted by , and 
the acceleiati(Hi vector, by f^ . The numbor of bade coirections is m. I'he coefficients |a^ 
and {0.J are determined in such a way that the highest possible order is obtained for a 
given k. It should be noted that for m = 0, the equations reduce to those of the Stormer- 
Cowell method. Therefore, the back-correction methods are simply a generalization of the 
classical Stormer-CoweU method. 

These methods can be used with pseudo-evaluations in the following algorithm, provided 
the acceleration can be separated according to dominant and perturbing terms: 

a. Predict a value, . for ^ i using equation 1 . 

b. Evaluate the dominant and perturbing accelerations using this value for ^ , saving 
the perturbing acceleration for subsequent calculations. 

c. Set 8 = 0. 

d. Using equation 2, obtain a corrected value, x® * ’ j , for , g . 

e. Reevaluate only the dominant acceleration using x® * ‘ .g and obtain f„ + , .g by 
adding the previously calculated perturbing acceleration at t^^ , g. 

f. If 8 = m, proceed to the next step of the integration. Otherwise, set 8 = 6 + 1 and 
go to (d) above. 

In problems involving earth satellites, the forces can readily be separated according to dom- 
inant and pertiubing terms. In such problems, these algorithms are especially efficient, 
because all evaluations of the forces after the first are simply pseudo-evaluations requiring 
only the reevaluation of the dominant forces (in this case, the two-body force). 

The methods using back corrections have been tested on several problems. A numerical 
int^ation of the orbit of the Applications Technology Satellite-F (ATS-F) has been per- 
formed in which case the errors were smaller by two or three orders of magnitude when 
compared to the classical Stormer-Cowell method using pseudo-evaluations. 

The results shown in table 1 are obtained when the seventeenth-order back -correction 
algorithm is applied to the to the Geodetic Earth Orbiting Sate!Iite-C (GEOS-C) orbit. 


( 1 ) 


( 2 ) 
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Also shown are the results obtained using the twelfth-orda' Stormer-Cowell algorithm 
presently avaflable in the Goddard Trajectory Determination System (GTDS). The errors 
are obtained by comparing with the results obtained using a step size of 40 seconds. It 
should be noted that for reasonable step sizes, the seventeenth-order back-correction method 
provides greater accuracy than the twelfth-order Stormer-Cowell method. Since pseudo- 
evaluations are used in both methods, the computer time required per step is the same 
(within a few percent) fOT both methods. 


Table 1 

Numerical Results for GEOS-C Satellite 


Step Size (s) 

Errors in Position after 24 Hours (km) 

\2th Order 
Stormer-Cowell 

nth Order 

Back-Correcticm Method 

60 

5.7 X 10 ® 

1.9 X lO-'' 

80 

3.2 X \(T* 

8.8 X 10^“' 

100 

1.1 X 10^^ 

1.5 X 10-^ 

200 

6.9 X 10-' 

2.2 X 10-* 


The stability region of the twelfth-order Stormer-Cowell method is, in fact, slightly smaller 
than that of the seventeenth-order back-correction method. A higher order Stormer-Cowell 
method would no doubt have exhibited greater accuracy than the twelfth-order Stormer- 
Cowell method, but the stability region would have been considerably reduced. Consequent- 
ly, the largest meaningful step size attainable with such a higher order method would have 
been much smaller. 

The methods using back corrections appear to be more efficient than the classical methods 
for problems in which the dominant and perturbing forces can be readily separated and in 
which the evaluation of the perturoing terms requires much more computer time than the 
evaluation of the dominant terms. 
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